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Electron  Current,  Beta  Limit  Line  Operation  and  Power  Balance  in  WB  Mode 


Introduction  and  Summary 

Analysis  of  electron  flow  in  the  outer  regions  of  the  Polywelltm  system  shows  the 
current  required  to  achieve  electron  /?  =  1  at  r  =  R  (the  <rfa>  =  1  condition).  This  is 
found  to  depend  on  the  square  root  of  the  injection  energy,  the  <rfa>  =  1  power 

required  will  vary  as  (Eo)1,5.  With  this  as  the  drive  power  an  expression  for  gross  power 
gain,  can  be  obtained  from  the  ratio  of  fusion  power  (integrated  over  system  volume) 
to  drive  power.  is  calculated  for  several  fuels.  Each  fuel  combination  has  unique 
bremmstrahlung  radiative  characteristics.  The  inclusion  of  bremmstrahlung  losses 
decreases  the  overall  gain,  GQa,  that  can  be  attained  by  any  system  to  values  less  than  Ggr, 
and  the  maximum  possible  value  for  G^  is  found  to  be  Ggr/2.  The  overall  gain  is  derived 
and  numerical  and  graphical  examples  are  given  to  show  the  parametric  range  of  system 
performance  with  such  losses. 


Results  show  that  DT  is  superior  and  easily  able  to  produce  high  gain  in  small 
systems  (e.g.  R  <  1.0  m).  It  is  also  found  that  DD  and  D  He  tend  to  give  comparable 
performance  for  comparable  mixture  ratios  (0.25  —  0.5),  and  that  both  are  considerably 
worse  (ca.  10  —  100x)*  than  DT  and  better  (10  -  lOOx)  than  p11B.  However,  D3He  is  seen 
to  be  no  better  than  pnB  if  operated  at  mixture  conditions  that  make  it  as  non— hazardous 
(radiation-free)  as  p^B.  These  analyses  suggest  that  there  is  little  incentive  for  use  of 
D3He  in  Polywelltm  systems.  They  also  show  that  pnB  CAN  be  made  to  operate  at  net 
power,  if  core  electron  energies  can  be  kept  sufficiently  low  (small  anode  height)  in  contrast 
to  conventional  magnetic  confinement  systems,  in  which  p1JB  can  NOT  be  made  to  yield 
net  power  (vs.  bremmstrahlung)  under  any  conditions. 

4 

*  TV'*  measure  of  "goodness"  used  here  is  the  functional  =  (BqR)  ft^  (see  later  text). 
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Electron  Current  Flow  at  Surface 0—1  Conditions 

Consider  a  system  operating  such  that  electron  0  =  1  at  all  times  at  the  edge  of  the 

n 

device,  thus  <*b>  =  (rb/R)  =  1.  At  this  condition  the  parameters  Z  =  8imcrc  and 
W  ss  Eo/(BqR)2  are  related  by  (ZW)  =  l.1  To  calculate  the  electron  injection  current 
required  to  drive  the  device  to  and  at  this  condition  it  is  neces^ry  first  to  know  the 
electron  density  njR)  =  nR  at  the  system  edge  r  =  R. 

The  electron  and  ion  charge  densities  are  essentially  equal  throughout  the  body  of 

1* 

the  device,  out  to  that  radius2  rk  at  which  the  surface  (edge)  "sheath"  begins.  At  this 
position  the  ion  density  and  electron  density  both  increase  with  increasing  r  -♦  R,  but  the 
electron  density  increases  more  rapidly  than  the  ion  density.  From  the  core  radius  r  =  r 

A 

to  r  =  r,  the  density  of  both  species  varies  closely  as  the  inverse  square  n(r)  =  n  (r  /r)  so 

K  C  C 

that  the  density  of  electrons  at  r  =  rfc  is  given  (approximately,  see  ref.  2)  by 

*(rk)  =  \  =  nc(rc/rk)2- 


The  total  radial  electron  current  at  rk  is 

’.,k  =  4mkVk  (1) 

where  vk  is  the  electron  radial  speed  at  rk>  The  radial  current  at  the  surface  at  r  —  R  must 
be  less  than  this,  because  conservation  of  magnetic  moment  and  associated  transverse 
momentum  in  the  electron  flow  forces  an  increase  in  transverse  electron  kinetic  energy  with 
increasing  radius  beyond  rk  (as  well  as  an  increase  inside  of  rk  to  the  "stagnation  radius" 

i* 

This  radial  position  ia  that  at  which  the  kinetic  energies  of  radial  motion  of  ions  and 
electrons  are  approximately  equal  (see  e.g.  ref.  2) 
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rf),  and  thus  a  decrease  in  radial  electron  energy  and  hence  in  radial  current  at  the 
boundary.  This  is  accounted  for  here  by  a  functional  factor  Fe(f1,<rk>)  <  1  that  depends 
on  the  degree  of  transverse  flow  of  the  electrons.  Thus  the  two  currents  can  be  related  by 


J«R  =  W.  =  4ltRVn  =  4ltIk VkF. 


(2) 


Now,  for  operation  at  <rfa>  =  1,  there  is  no  mirror  reflection  effect  and  the  electron 
energy  will  simply  follow  the  shape  of  the  internal  potential  well  which,  in  turn,  will  follow 
that  of  the  internal  B  field.  This  will  follow  the  functional  form  of  the  original  unperturbed 
field,  even  though  the  field  is  at  reduced  amplitude  due  to  the  diamagnetic  current 
reduction  of  field  strength  by  the  factor  fQ  =  B^JBq(t)  within  the  wiffle  ball  (WB) 
surface.3  With  this  shape,  it  can  be  shown  that  the  total  electron  kinetic  energy  follows 
the  formula1 


Ek(<r>)  =  Eo<r>mfo(<r>)  (3) 

for  a  rollover  field  of  index  exponent  (m).  For  a  truncated  cube  m  =  3,  the  rollover 
function  is  fQ  =  2/(l+<r>5),  and  thus  the  velocities  vary  about  as 
T«rk»  =  v(l)[2<rk>J/(l+<rl>5)]0-5. 

The  functional  factor  F#  has  been  derived  previously4  in  terms  of  the  "stagnation 
radius"  rf  for  electron  flow  in  the  outer  regions  of  the  system,  where  electrons  establish  the 
potential  well  that  accelerates  the  ions.  Adjusting  this  to  the  "critical  radius"  rk  by  a 
correction  term  for  the  difference  in  potential  and  throughflow  area  between  the  stagnation 
radius  <rp  =  (dEj^E^)0-5  =  (fj.)0'5  and  the  critical  radius  <rfc>  at  which  the  ion  and 
electron  density  gradients  change  sign,  and  taking  account  of  diamagnetic  field  reduction 
within  the  WB  surface,  gives  this  as  approximately 
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Note  that  the  diamagnetic  factor  fD  cancels  out  of  the  numerator  of  eq.  (4);  it  is  shown 
simply  for  completeness  of  description.  With  this,  and  including  the  multiplying  effect  of 
electron  recirculation  by  the  factor  G^,  the  surface  electron  injection  current  is  then  simply 


(5) 


for  <rfa>  =  1  (r  =  R),  where  k  =  6.28E18  chgs/Asec  converts  to  current  Ie  in  amps,  the 
injection  speed2*  of  electrons  is  v„  =  (2E  /m  )*,  and  all  units  are  cgs. 

K  O  C 


Note  that  this  speed  is  that  at  which  the  electrons  are  initially  injected  purely  radially,  while 
the  actual  electron  motion  at  and  near  the  boundary  of  the  device,  after  the  first  pass  through  the 
system,  will  be  very  tangential  due  to  the  conservation  of  transverse  momentum  and  of  magnetic  moment 
in  electron  flow  through  the  multicusp  system,  as  mentioned  earlier.  The  functional  factor  corrects 
for  this  use  of  a  "false"  radial  speed  at  r  =  R.  Note  that  this  situation  holds  true  for  0  =  1  operation 
out  to  any  radius,  as  beyond  or  directly  "at"  the  0  =  1  radius  the  external  B  field  is  still  imposed  at  its 
original  strength,  and  the  electrons  will  still  "mirror"  and  conserve  momentum  through  the 
magnetic  moment  effect  in  these  regions,  as  well  as  inside  the  WB  sphere  if  enough  residua]  field  is  left 
therein.  Finally,  it  is  important  to  note  that  operation  at  0  =  1  will  almost  certainly  be  cut  off  at 
r  =  r^,  due  to  the  0  instability  beyond  this  radius,  thus  the  region  between  r^  and  R  will  still  act  in  an 
MR  mode  fashion. 


4 


Electron  Injection  Power  on  the  <rfa>  =  1  line 

Taking  (ZW)  =  1  on  the  <rfa>  =  1  line  and  noting  that  =  2reZ/N(kLS)2(l-aR) 
2 

here  and  Z  =  8xncrc,  allows  reduction  of  eq.  (5)  to  yield 


(fo<ri.>\)0'5  ”(''rlS)i(i-°fe) 
(2mt)°-5(2tekt) 


(PKV 


0.5 


(6) 


Here  r  =  2.828E— 13  cm  is  the  electron  Compton  radius,  N  is  the  effective  cusp  number,  k 
c  L 

is  the  effective  loss  radius  factor^  (on  surface  gyro  radius),  S  is  the  average  sine  of  the 
electron  collision  angle  with  the  remaining  cusp  field  at  the  adiabaticity  radius, 

=  agNg/N  is  the  effective  repeller  effectiveness  ay  for  Ng  guns  each  with  af  =  agl  and 
ke  =  1.6E-12  ergs/eV  allows  Eq  to  be  used  in  units  of  eV. 


As  an  example,  consider  a  truncated  cube  system  with  typical  values  of  f±  =  0.3 

(wide  divergence  of  electrons),  <rk>  =  0.83  and  <rc>  =  IE— 2.  With  these,  the  functional 

factor  to  correct  for  inertial  flow  convergeqce  and  magnetic  and  transverse  momentum 

conservation  effects  becomes  Fe  =  0.0610.  With  this,  and  taking  N  =  8.0,  Ng  =  6  with 

a  =  0.9,  k  =  2.0  and  S2  -  0.6  (mean  of  range  from  ir/6  to  2/3), 5  the  electron  drive 
©  u 

current  required  is  found  to  be  related  to  the  injection  energy  by 

IeR  =  2.88(E/5  (7) 

for  operation  along  the  <rfe>  =  1  line.  Here  Eq  is  in  eV. 


Thus  if  Eq  =  1E4  eV,  for  example,  the  drive  current  must  be  IeR  =  288  A;  if 
Eo  =  2.5E4  eV  (25  keV),  then  IeR  =  455  A. 


5 


I 


The  <rb>  =  1  line  can  be  reached  most  directly  by  approaching  WB  mode 

operation  from  a  small  value  of  G.  (e.g.  along  the  line  of  G.  ,  =  1),  starting  operation  at 

J  jw® 

small  current  and  increasing  current  until  <rb>  =  1  is  attained.  This  requires  that 


G.  .E  2r 

g.  ,w=  >*b  ;= 


jwb 


(B.H)' 


N(k,S)‘ 


(8) 


Using  this  in  eq.  (6)  shows  that  the  critical  injection  current  at  G^  =  1  is  related  to  the 
system  field  strength  and  size  by 


*eR  ” 


N 


r  m 
e  e 


f  <r,  > 

o  k 


0.5 


2k 


F(B0R) 


(9) 


For  the  previous  parameters  this  becomes  IeR  =  0.39(BoR)  A,  so  that  the  field,  size  and 
injection  energy  are  related  by  (BqR)  =  7.37(Eq)0'5.  For  R  =  92  cm,  for  example,  this 
yields  Bq  =  8.01  G  for  E  =  1E4  eV,  IeR  =  288  A  as  illustrated  above,  or  directly  from 
eq.  (8).  At  this  condition  the  parameter  W  =  Vfl  —  0.9064E— 13  cm,  or 
WA  =  56.65  keV/(kGcm)2. 


For  Eq  =  1E4  eV,  as  before,  PeR  =  2.88  MW  is  required  along  the  <rb>  =  1  line. 
Having  reached  this  line  at  G^wb  =  1  it  is  then  possible  to  move  up  along  this  line  by 
reducing  the  value  of  W  from  that  at  Wj  to  successively  lower  values,  keeping  the  injection 
current  constant  (or  very  nearly  so).  On  this  line  G^wb  will  vary  inversely  with  W;  thus  * 
raising  the  B  field  by  lOx  will  increase  G^wb  by  lOOx,  and  increase  Z  by  lOOx. 

Finally,  note  that  this  point  (GjWB  =  1»  <rb>  =  *)  the  <rb>  =  1  ^ne  can 
attained  only  by  injection  of  the  required  current  at  an  electron  drive  power  of 
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PeR  =  2.88(Eo)1,5  watts.  Any  system  to  produce  net  power  with  these  parameters  must 
operate  with  more  than  2.88  MW  of  fusion  power  generation. 

Fusion  Power  Generation  and  Gross  Gain 


The  fusion  power  output  is  given  as  the  core-generated  power  multiplied  by  a 
correction  factor,  in  the  approximate  algorithmic  form 


Pfus  - 


4m 


(Kf)(nfafo[Eo(  1-r7)]i)(Ef)-f- 


k1*5 
\  e 

•2Eo(l-,)- 

) 

10 

m. 

L  i 

0.5 

b.. 

u 


(10) 


Here  the  factor  Kf  =  2(1 — 77q) /( 1— ry)  accounts  for  the  fact  that  up  to  0.75-0.85  of  the  power 
is  generated  outside  the  convergence— limited  core  at  radius  r  ,  depending  upon  the  anode 
height,  7].  ®  Here  the  fusion  cross-section  has  been  taken,  for  illustrative  purposes,  to  be 
ufu#  =  <rfoE®  over  the  energy  range  of  interest.4  This  simple  formulation  allows  analytical 
illustration  of  the  relationships  between  fuel  properties,  system  parameters,  and  system 
gross  gain.  More  exact  cross-section  formulae  are  used  later  for  computer  calculations  of 
these  relationships.  Recalling  that  nc  =  n^  +  n2Z2,  where  nt  =  fjn.  and  n2  =  f2n.,7  the_# 
core  density  is  found  from  the  identities  ZW  =  1,  Z  =  8imcrc  and  W  from  eq.  (8),  to  be 


n.  = 

i 


N(kLS)2Gi.b 

2r  8 7rr  2 
e  c 


1 

i+yvij 


(11) 


With  this  the  fusion  power  becomes 


Pfu.= 


bu  V^a/S)11  ■ 8  ( £.(1-1,)]'+  °  '5(tLS)4(V.)15(N2) 


3  x  64  x  10t  (re)\[l+f2(Z2-l)] 


(12) 
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Since  G^b  =  Wj/W  =  (B^/Bj)2  for  Wj  taken  at  Gjwb  =  1,  <rb>  =  1.  it  is  evident 
that  Pfu>  scales  as  B*.  Note  that  this  is  the  same  as  and  directly  analogous  to  the  power 
scaling  of  all  other  magnetic  confinement  fusion  concepts.  The  power  amplification  or 
gross  power  gain  G^  is  this  expression  eq.  (12)  divided  by  IeREo  from  eq.  (6),  or 

G  W.  ( me/5.)°»(KfgfoEtN)(k1S)2(Eo)-1^) 

,r  (fo<rk>3)0'5(480nrc)Ft[l+f2(Z2-l)];!(l-aR)  jwb 

a 

for  Ef  in  MeV  and  Eq  in  eV,  with  all  other  units  in  cgs,  but  afo  =  cm  /(eV)  . 

For  DT  fuel  Ef  =  17.6  MeV  and  b„  =  0.25  (for  50:50  mixtures).  At  Eq  =  2E4  eV 
(20  keV),  the  fusion  reaction  cross-section  varies  about  as  s  =  2.7  with  a  coefficient  given 
approximately  as  afo  =  3.5E-36  cm2/(eV)2'7.  Taking  rj  =  0.1,  m.  (reduced)  =  1.2  mp  and 
r.  as  1.0  cm,  and  other  parameters  as  before,  eq.  (13)  gives 

Gg,lDT  =  2-60E-10(Gjwb)2  (14) 

Thus  net  power  can  be  made  only  if  G^wb  >  6.2E4.  For  pnB  (50:50,  with  in.  =  0.9  mp)  at 
300  keV  (s  =  2.85,  afo  *  4.48E-41  cm2/(eV)2,85)  it  is  similarly  found  (for  rc  =  1.0  cm) 
that 


^gJpllB 


1.76E— 13(Gjwfa)2 


(15) 


so  that  net  power  requires  that  G^wb  >  2.4E6.  For  a  system  gross  power  gain  of  G^  =  40 
with  DT,  G.wb  >  3.93E5,  while  for  pnB,  G^b  >  1.51E7,  for  the  parameters  chosen  in  the 
foregoing  examples. 
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It  is  instructive  to  restate  these  gross  gains  in  terms  of  the  system  size  and  fields 
required.  To  do  this,  return  to  eq.  (8)  and  write  G^b  in  terms  of  (BqR)  and  Eq1  and 
substitute  into  eqs.  (13,14,15),  thus 


k  r 

•  e 


(mt/5.)08(b..Kf.rtt|Et)(Eo)-"(l-^) 
(fo<rk>  3  )°  •  5(120«c)  (  N  )(kLS)2Ft[l+f2(Z2-l ) 


■-3/ 


it+0.5 


(V) 


(16) 


and 


Ggr  I DT  =  0022E-20(BoR)4/tc  (17a) 

for  Bq  in  gauss  (G),  R  and  rc  in  cm,  or 

G„lDT  =  a022<BoR)X  (17b) 

for  B  in  kG,  R  and  r  in  m.  For  example,  for  <r  >  =  0.01,  if  R  =  0.5  m,  B  =  4.5  kG, 

o  c  c 

then  G  =  1.13;  if  R  =  0.8  m,  B  =  8.0  kG,  then  G  =  46.1.  And  for  pnB  at  the  previous 
gr  gr 

conditions 


g,lpllB  =  0066E-26(BoR),/'c  (C«5) 

(18a) 

«,lpilB  =  0  °66E-6(b0r)4Ac  (K3.“) 

(18b) 

Thus,  for  <r  >  =  0.01,  if  B  =15  kG  then  R  =  8.0  m  is  needed  for  G  =  1.14.  If 

*  c  ’  o  gr 

Bq  =  25  kG  and  R  =  10.0  m,  then  the  gross  gain  can  be  Ggr  =  17.2. 
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Note  from  eq.  (16)  that  will  increase  with  increasing  Eq  only  if  the  fusion 

A 

cross-section  changes  more  rapidly  than  s  =  3.  Since  p  Li  and  DD  reactions  (as  well  as 
some  others)  have  s  <  3  over  most  of  their  energy  range  of  interest,  it  is  dear  that  high 
energy  operation  is  not  necessarily  optimal  for  use  of  such  fuels. 

Bremmstrahlung  Losses  and  System  Overall  Gain  in  <rfa>  =  Of  ode 


However,  the  above  gain  formula  (eq.  16)  does  not  account  for  bremmstrahlung 
losses.  If  these  are  induded  it  is  found  that  low  energ-  operation  leads  to  excessive 
bremmstrahlung  vs.  fusion  power  density.  There  is  thus  an  optimum  value  for  well  depth 
(mean  reaction  energy)  for  any  fuel  combination  with  s  <  3.  In  the  region  where 
cross-section  variation  is  8  >  3,  there  is  no  optimum  well  depth  driven  by  bremmstrahlung, 
and  the  best  well  depth  is  that  at  the  highest  possible  cross-section  for  s  >  3,  so  long  as 
the  core  electrons  remain  cold.  Calculating  bremmstrahlung  power  as  core-generated 
radiative  power  multiplied  by  a  correction  factor  Kb,  for  extra-core  bremmstrahlung  (as 
for  fusion  power),  the  fusion/bremmstrahlung  power  ratio  is  given  from  eq.  (10)  and  the 
standard  formulae  for  core  bremmstrahlung  ^mission,7  evaluated  at  core  conditions,  as 


<pn>>  = 


,-0'5(2/m.)0,5b..k  1S 
_ v  '  V  ij  e 


1.69e-31  K.F  rj 
b  > '« 


0.5 


(19) 


where  =  [l+f2(Z2— 1 )] [  1 +f2( Zj — 1 )] •  The  factor  Kb  is  found  to  be  approximately  2.5  for 
most  systems  of  interest.  With  this,  the  system  overall  gain  (G^),  induding 
bremmstrahlung  losses  (but  not  induding  ohmic  power  to  drive  the  magnetic  field  coils), 
will  be 
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1 


1 


(20) 


G  = 

oa 


G  <P,k> 

gr  fb 


Writing  Ggf  =  Ce(Eo),_3  (from  eq.  16)  and  <Pft>  =  Cb(Eo)*  (from  eq.  19)  and 
differentiating  eq.  (20)  to  find  the  optimum  condition  for  Eq  for  those  fuels  with  s  <  3 
yields 


o' opt 


sC 


(3-s)C. 


1/3 


(21) 


and  if  s  =  1.50,  for  example  (a  reasonable  fit  to  several  of  the  fusion  reactions  over  the 
midrange  of  their  ascending  cross-section  energy  region),  then  EJ  =  and 

the  overall  gain  eq.  (20)  gives  the  maximum  possible  gain  as 


g  i 

oa'  max 


c  c, 

e  b 


(22) 


For  s  =  1.50,  eqs.  (20)  and  (21)  require  that  G  =  <P_>  at  optimum  operation. 
At  this  condition  eq.  (22)  shows  that  GoJmax  =  Ggr/2  is  the  largest  possible  overall  gain 
of  the  system.  Thus,  the  maximum  performance  of  a  fuel  with  s  =  1.5  will  be  set  by  the 
value  of  Frr  =  (BQR)4/rc  =  B4R3/<rc>;  larger  values  of  (BqR)  give  higher  values  of  Goa- 
All  cross-sections  eventually  reach  a  peak  with  increasing  energy,  at  which  point  s  =  0. 
Thus  as  the  reaction  energy  (well  depth)  is  increased  above  the  region  where  s  =  1.5,  this 
exponent  decreases  steadily  to  8  =  0,  and  becomes  negative  thereafter.  At  energies  abc  *e 
the  s  =  1.5  region,  the  af/E^  term  decreases  more  and  more  rapidly  and  the  effect  of 
increasing  FRR  becomes  less  and  less  until,  as  <7f  approaches  its  peak  value,  increasing  FBR 
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has  no  effect  at  all.  This  behavior  is  seen  quite  clearly  in  graphic  portrayals  of  these 
parameter  variations. 

These  can  be  obtained  from  the  foregoing  analyses  in  considerably  improved  form  by 
use  of  correct  functional  formulae  for  the  fusion  cross-6ections  of  the  various  candidate 
fuels.  For  convenience  these  are  summarized  in  Table  1,  for  several  fuels.  The  table  gives 
the  cross-section  variation  with  mean  collision  energy  in  the  (^t/CG  system  that 
characterizes  these  Polywelltm  converging-flow  systems.  Using  the  functional  form  cr^E), 
where  E  =  Eo(l-j;),  and  letting  =  0,  recasts  eq.  (16)  as 

G  _  (V  «  A  ,  )[(l-?)me/j5.) 0  -5(  b.  .K  f  Ef )  [  gf(E)l _ 

V  (fo<rk>3)0S(120T<rc>)(Eo)3(N)(kLS)2Fe[l+f2(Z2-l)]2  ° 

Figures  l(a-d)  show  the  parameter  Ggr/FBR>  where  FRR  =  (BoR)4/rc,  as  a  function 
of  electron  injection  energy  Eo,  with  fx  =  0.3,  for  the  fuels  in  Table  1,  for  several  values  of 
the  anode  height  parameter  tj  and  a  core  convergence  ratio  of  <t^>  =  0.01,  for  equal 
fractional  mixtures  of  fuel  components,  f,  =  f2  =  0.5  +  f2  =  1-0),  where  fj  is  the 

fractional  mix  of  singly-ionized  (Z  =  1)  fuel  and  f2  is  that  of  the  multiply-ionized  (Z  >  1) 
fuel  component. 

Figures  l(e-h)  show  Ggr/FfiR  for  fuel  fraction  f2  at  optimum  values  to  yield  the 
maximum  ratio  of  fusion  to  bremmstrahlung  power,  taken  from  prior  study  of 
bremmstrahlung  losses.7  These  are  f2  =  0.261  for  D3He,  0.084  for  pnB  and  0.5  for  DT.  • 
The  curves  for  DT  and  DD  are,  of  course,  identical  in  both  sets  of  figures.  Note  that  DD 
surpasses  D  He  in  gross  gain  performance,  that  DT  is  clearly  best,  and  p  B  most  difficult 
and  least  promising. 
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*  2  * 

The  value  of  G  scales  simply  with  <r  >  as  G  =  G  (10  /<r  >)  where  G  is 

gr  r  *  c  gr  gr'  '  c  '  gr 

the  gross  gain  value  as  shown  in  Figures  1.  It  is  interesting  to  note  the  position  of  the 
maxima  for  the  different  fuels.  Increasing  anode  height  results  in  lower  values  of  gross 
gain,  and  shifts  the  peak  to  slightly  higher  Eq  values.  Also  note  that  the  gross  gain  for 
D  He  and  p  B  is  only  slightly  improved  at  optimum  bremmstrahlung  mixture  ratio  than 
for  50:50  mixtures. 

Use  of  this  equation  together  with  eq.  (19)  for  bremmstrahlung  allows  exact 
computation  of  the  overall  gain  including  bremmstrahlung  G  from  eq.  (20).  This  is 
shown  in  Figures  2-5  as  a  function  of  Eq  for  each  fuel,  for  a  range  of  the  parameter  FfiR, 
and  for  two  values  of  the  minimum  core  electron  energy  (anode  height),  rj  =  77,  of  0.01  and 
0.1.  The  bremmstrahlung  extra-core  factor  was  taken  to  be  =  2.5  in  these  calculations. 

Figures  2(a,b)  and  3(a,b)  show  overall  gain  for  DT  and  DD,  respectively,  for 
f l  =  f2  =  0.5  mixtures  at  two  values  of  j?e  =  rj.  Note  that  the  position  of  maximum  GQa 
has  shifted  upward  in  energy  very  slightly  for  DT,  but  up  by  as  much  as  lOx  for  DD;  a 
result  of  the  weak  energy  dependence  of  the  DD  fusion  cross-section  to  counter  balance 
bremmstrahlung  losses. 

»  W 

Figures  4(a-d)  and  5(a-d)  show  G^  for  D3He  and  pnB,  respectively,  for  50:50 
mixtures  (Figs.  a,b)  and  bremmstrahlung-optimum-f2  mixtures  (Figs.  c,d)  at  two  values 
of  rj'  =  rj.  For  both  fuels  the  gain  is  increased  and  the  peak  injection  energy  required  is 
higher  at  the  optimum  f2  value.  The  effect  with  p  B  is  more  pronounced  than  for  D  He, 
and  it  is  clear  that  pnB  can  give  net  power  only  at  small  values  of  rj  <  0.03,  no  matter 
how  large  FRR  may  be  made. 
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As  an  example  of  the  situation  consider  a  pnB  system  operating  at  Eq  =  500  keV 
with  G  =  7  9  at  F  =  1E11  =  (B  R)4/r  .  If  <r  >  =  IE-3,  for  example,  the 
B4R3  =  1E8  (kG)4(m3).  Thus,  if  Bq  =  31.6  kG,  R  =  4.64  m  is  required,  or  if  R  =  2.5  m, 
then  Bq  =  50.3  kG  must  be  supplied.  These  are  large  parameters  by  the  standards  of  the 
small  systems  attainable  with  DT,  but  still  small  compared  with  magnetic  confinement 
fusion  systems.  The  fact  that  pnB  can  be  made  to  yield  ne^ower  is  a  result  of  the  unique 
energy  and  particle  distributions  inherent  in  the  Polywelltm  inertial-electrostatic  system. 
This  useful  result  is  in  striking  contrast  to  all  conventional  magnetic  confinement  systems, 
where  pnB  can  not  be  made  to  work  at  all. 

4 

Comparison  of  DD  and  D  He  shows  that  the  two  offer  comparable  performance 
(Gq4)  at  Fbr  *s  1E8  but  that  the  overall  gains  for  D3He  become  larger  than  those  for  DD 
at  Fbr  >  1E8;  a  factor  of  3x  higher  at  FfiR  =  1E9.  However,  the  injection  energy  (well 
depth)  is  about  2x  higher  at  this  latter  condition. 

Finally,  while  it  is  clear  that  D  He  offers  much  higher  performance  than  p  B  for 
the  mixtures  assumed,  at  smaller  values  of  both  E  and  Fnt>,  it  is  important  to  note  that 

O  DA. 

this  performance  includes  a  very  large  neutron  production  rate  and  consequent  radiation 
hazard,  not  present  with  pnB. 

Figure  6  shows  this  neutron  production  rate  in  terms  of  the  DD  contribution  to  total 
fusion  power  in  a  fj  =  f2  =  0.5  (50:50)  D  He  mixture  at  various  well  depths.  One-half  of 
the  DD  reactions  produce  a  2.45  MeV  neutron.  Note  that  the  DD  power  fraction  never 
drops  below  about  0.025  of  the  total,  and  that  operation  at  well  depth  above  this  minimum 
will  yield  increased  neutron  power.  Of  course,  in  the  Polywelltm  system  this  can  be 
reduced  by  changing  the  fuel  mixture  ratio,  increasing  He  and  decreasing  D. 
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The  fusion  power  goes  roughly  as  nDnHe  while  the  DD  (neutron)  power  goes  as 
(nD)2,  thus  the  fractional  DD  neutron  power  will  vary  as  nHe/nD  =  ^/(l— f2)-  Reduction 
of  neutron  radiation  then  requires  that  f2  be  made  large  (i.e.  close  to  1.0).  If  f2  =  0.999,  for 
example,  the  neutron  power  will  be  down  by  about  1E6  from  that  of  the  50:50  mixture 
condition  and  the  system  will  not  require  any  significant  radiation  shielding.  At  this  level, 
D  He  is  comparable,  in  terms  of  radiation  hazard,  to  p  B. 

However,  this  mixture  is  far  off  of  optimum  for  bremmstrahlung  losses,  which  will 
be  at  nearly  their  maximum  possible  level  -  such  a  system  is  essentially  "all"  He  with 
Z  =  2.  Applying  the  previous  formalism  to  this  (1:1000)  mixture  (f2  =  0.999)  gives  the 
performance  curves  shown  in  Figure  7.  From  this  it  is  clear  that  "neutron-free"  D  He 
systems  offer  less  performance  potential  than  does  pnB  and,  in  fact,  they  can  not  yield  net 
power  at  all  with  any  attainable  values  of  FRR. 

9 

The  utility  and  potential  value  of  He  as  a  fusion  fuel  thus  may  be  open  to  question. 
D  competes  with  it  as  a  neutron-producing  fuel,  and  nB  out-performs  it  in  a 

9 

"neutron-free"  mode.  He  is  very  scarce  and  will  remain  costly  even  if  mined  from  the 
lunar  surface,  while  nB  is  readily  available  at  reasonable  cost  on  the  Earth,  and  D  is  found 
in  practically  unlimited  quantities  in  all  water  and  is  virtually  free.  However,  intermediate 

9 

D  He  fuel  mixtures  (e.g.  1:30)  allow  for  a  substantial  reduction  in  shielding  over  DD 
systems  due  to  the  lower  neutron  rate  and  run  at  lower  injection  energies  than  comparable 
p  B  systems.  Therefore,  the  use  of  D  He  in  this  intermediate  range  needs  to  be  further 
explored  to  determine  its  true  value. 

Of  course,  the  power  gain  scalings  given  above  do  not  include  the  power  required  to 
supply  the  Bq  fields,  thus  they  are  reasonably  correct  only  for  superconducting  magnet 
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systems.3*  If  magnets  are  of  normal  conductors  (e.g.  copper  and/or  its  alloys)  the  systems 
will  be  driven  to  still  higher  injection  energies  and  larger  sizes  for  optimum  operation. 


Finally,  analyses  of  the  systemic  effects  of  synchrotron  losses  have  also  been  made 
and  are  reported  elsewhere.®  These  show  optimum  operating  conditions  (i.e.  mixture 
fractions  and  EJ  )  different  than  those  for  bremmstrahlung  optimization.  Since 
synchrotron  radiation  will  be  partially  self-absorbed  and  can  be  partially  reflected,  greater 
flexibility  in  system  design  is  allowed  for  its  control.  For  this  reason  (i.e.  it  can  be  made 
"design-specific’')  it  has  not  been  included  here.  Bremmstrahlung  losses  can  not  be 
mitigated  by  system  design  changes,  they  are  endemic  to  the  plasma  mixture,  thus  are 
included  in  the  analyses  of  inherent  gain  given  here. 


Studies  of  the  effects  of  variable  D3He  mixtures  and  of  magnet  drive  power 
requirements  are  currently  underway  and  will  be  reported  in  future  EMC2  Technical 
Notes/Reports. 


But  sot  exactly,  because  s/c  magnets  require  cryogenic  cooling  power  for  their  operation  and 
this  power  must  be  included  in  the  power  balance  for  Goa-  However,  the  s/c  cooling  power  is  only 
weakly  coupled  to  the  B  field  strength  thus  the  effect  is  less  profound. 
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TABLE  1 
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Figure  1.  (*— d) 

4 

Gross  gain  in  terms  of  FgR1  where  FgR  =  (B^R)  /r^,  as  a  function  of  electron  injection  energy, 
E^,  for  the  fuels  listed  in  Table  1.  Fuel  ratios  are  50:50,  plotted  for  a  range  of  anode  height 
parameter,  0.01  <  TJ  <  0.3. 

Figure  1.  (e— h)  ^ 

Gross  gain  in  terms  of  Fodi  where  F__  =  (B  R)  /r  ,  as  a  function  of  electron  injection  energy, 

Da  da  O  C 

Eq,  for  the  fuels  listed  in  Table  1.  Fuel  ratios  are  for  optimum  <P^>,  plotted  for  a  range  of 
anode  height  parameter,  0.01  <  Tf  <  0.3. 

Figure  2.  (a,b) 

System  overall  gain  for  DT  for  anode  height  parameter  T)  —  0.01  and  0.1,  plotted  as  a  function  of 

4 

electron  injection  energy  for  a  range  of  F__  =  (B  R)  /r  ,  1E3  <  F__  <  1E6. 
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Figure  3.  (a,b) 

System  overall  gain  for  DD  for  anode  height  parameter  T)  =  0.01  and  0.1,  plotted  as  a  function  of 

4 

electron  injection  energy  for  a  range  of  FflR  =  (BqR)  /r^  1E6  <  FfiR  <  1E9. 

Figure  4.  (a-d) 

3 

System  overall  gain  for  D  He  for  anode  height  parameter  7]  =  0.01  and  0.1,  plotted  as  a  function 

4 

of  electron  injection  energy  for  a  range  of  FRR  =  (BqR)  /r^,  1E6  <  FRR  <  1E9  for  both  50:50 
(a,b)  and  bremmstrahlung— optimised  (c,d)  fuel  mixtures. 

Figure  5.  (a— d) 

System  overall  gain  for  p  B  for  anode  height  parameter  T)  =  0.01  and  0.1,  plotted  as  a  function 

4 

of  electron  injection  energy  for  a  range  of  P--  =  (B  R)  /r  ,  1E8  <  F_„  <  1E11  for  both  50:50 

jdK  o  c  rJK 

(a,b)  and  bremmstrahlung— optimised  (c,d)  fuel  mixtures. 

Figure  6. 

3 

Percentage  of  total  D  He  fusion  power  due  to  DD  side  reactions  as  a  function  of  electron 
injection  energy  for  50:50  fuel  mixture.  Plotted  for  the  range  of  anode  height  parameter, 

0.01  <  TJ  <  0.3,  used  in  Figures  1. 

Figure  7. 

3  3 

"Neutron— free"  D  He  system  performance,  with  fj  ss  0.999  (D:  He  =  1:1000). 


18 


REFERENCES 


1  Bussard,  R.W.,  K.E.  King,  "Electron  Recirculation  in  Electrostatic  Multicusp  Systems:  II  —  System 
Performance  Scaling  of  One— Dimensional  ’Rollover*  Wells",  Energy/Matter  Conversion  Corporation 
Report  EMC2— 0791—04,  July  1991 

3  Bussard,  R.W.,  "Ion  and  Electron  Flow  and  Some  Critical  Radii  in  PolywelltB  Systems", 
Energy/Matter  Conversion  Corporation  Report  EMC2— 0791— 02,  July  1991 

3  Bussard,  R.W.,  "Effective  Gyro  Hole  Loss  Radius  and  Diamagnetic  Lkrut  in  Polywelltm  Systems", 
Energy/Matter  Conversion  Corporation  Report  EMC2— 0591— 02,  May  1^1 

4  Bussard,  R.W.,  G.P.  Jellison,  G.E.  McClellan,  "Preliminary  Research  Studies  of  a  New  Method  for 
Control  of  Charged  Particle  Interactions",  Pacific— Sierra  Research  Corp.  Report  PSR  1899,  30  November 
1988,  Final  Report  under  Contract  DNA001— 87— C— 0052 

5  Bussard,  R.W.,  K.E.  King,  "Electron  Recirculation  in  Electrostatic  Multicusp  Systems:  I  — 
Confinement  and  Losses  in  Simple  Power  Law  Wells",  Energy/Matter  Conversion  Corporation  Report 
EMC2— 0491— 03,  April  1991 

•  Bussard,  R.W.,  K.E.  King,  "Virtual  Anode  Height  Variation  in  PolywelltB  Systems",  Energy/Matter 
Conversion  Corporation  Report  EMC2— 0991— 02,  September  1991 

7  Bussard,  R.W.,  K.E.  King,  "Bremmstrahlung  Radiation  Losses  in  PolywelltB  Systems",  Energy/Matter 
Conversion  Corporation  Report  EMC2— 0891— 04,  August  1991 

8  Bussard,  R.W.,  K.E.  King,  "Edge  Region  Distributions  and  Synchrotron  Radiation",  Energy/Matter 
Conversion  Corporation  Report  EMC2— 0991— 04,  September  1991 


19 


(e-ui>~o>i/i)  [t-(a°a)/(o-i)p3o 


(C-uit-OH/T) 


/t\  r*„fuoa\ /(->  id  i»«. 


Electron  Injection  Energy’ (keV)  >  Electron  Injection  Energy  (keV) 

Figure  1.  (a-d)  ;  Figure  1  c.  Figure  1  d. 

Gross  gain  in  terms  of  FfiR,  where  FRR  =  (BQR)4/rc,  as  a  function  of  electron 
injection  energy,  Eq,  for  the  fuels  listed  in  Table  1.  Fuel  ratios  are  50:50,  plotted  for 
a  range  of  anode  height  parameter,  0.01  <  r/  <  0.3. 
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Figure  2.  (a,b)  Electron  Injection  Energy  [keV]  Figure  2  b. 

System  overall  gain  for  DT  for  anode  height  parameter  tj  =  0.01  and  0.1,  plotted  as 
a  function  of  electron  injection  energy  for  a  range  of  FBR  =  (BqR)  /rc, 

1E3  <  <  1E6. 
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Figure  3.  (a,b) 

System  overall  gain  for  DD  for  anode  height  parameter  r?  =  0.01  and  0.1,  plotted  as 
a  function  of  electron  injection  energy  for  a  range  of  Fbr  =  (BQR)4/rc, 

1E6  <  Fbr  <  1E9.  11 
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system  Uveraii  Gain  ror  DU  vs.  Injection  Energy 
for  indicated  values  of  (BoR)~4/rc  [kG~4m~3] 
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Figure  3.  (a,b) 

System  overall  gain  for  DD  for  anode  height  parameter  ij  =  0.01  and  0.1,  plotted  as 
a  function  of  electron  injection  energy  for  a  range  of  =  (BQR)4/rc, 

1E6  <  F«„  <  1E9. 
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